Helical motors and formins synergize to compact chiral filopodial bundles: A theoretical perspective.

Chiral actin bundles have been shown to play an important role in cell dynamics, but our understanding of the molecular mechanisms which combine to generate chirality remains incomplete. To address this, we numerically simulate a crosslinked filopodial bundle under the actions of helical myosin motors and/or formins and examine the collective buckling and twisting of the actin bundle. We first show that a number of proposed mechanisms to buckle polymerizing actin bundles without motor activity fail under biologically-realistic parameters. We then demonstrate that a simplified model of myosin spinning action at the bundle base effectively "braids" the bundle, but cannot control compaction at the fiber tips. Finally, we show that formin-mediated polymerization and motor activity can act synergitically to compact filopodium bundles, as motor activity bends filaments into shapes that activate twist forces induced by formins. Stochastic fluctuations of actin polymerization rates and slower cross linking dynamics both increase buckling and decrease compaction. We discuss implications of our findings for mechanisms of cytoskeletal chirality.

Computing hydrodynamic interactions in confined doubly periodic geometries in linear time.

We develop a linearly scaling variant of the force coupling method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly periodic geometry with either a single bottom wall or two walls (slit channel) in the aperiodic direction. Our spectrally accurate Stokes solver uses the fast Fourier transform in the periodic xy plane and Chebyshev polynomials in the aperiodic z direction normal to the wall(s). We decompose the problem into two problems. The first is a doubly periodic subproblem in the presence of particles (source terms) with free-space boundary conditions in the z direction, which we solve by borrowing ideas from a recent method for rapid evaluation of electrostatic interactions in doubly periodic geometries [Maxian et al., J. Chem. Phys. 154, 204107 (2021)]. The second is a correction subproblem to impose the boundary conditions on the wall(s). Instead of the traditional Gaussian kernel, we use the exponential of a semicircle kernel to model the source terms (body force) due to the presence of particles and provide optimum values for the kernel parameters that ensure a given hydrodynamic radius with at least two digits of accuracy and rotational and translational invariance. The computation time of our solver, which is implemented in graphical processing units, scales linearly with the number of particles, and allows computations with about a million particles in less than a second for a sedimented layer of colloidal microrollers. We find that in a slit channel, a driven dense suspension of microrollers maintains the same two-layer structure as above a single wall, but moves at a substantially lower collective speed due to increased confinement.

Bending fluctuations in semiflexible, inextensible, slender filaments in Stokes flow: Toward a spectral discretization.

Semiflexible slender filaments are ubiquitous in nature and cell biology, including in the cytoskeleton, where reorganization of actin filaments allows the cell to move and divide. Most methods for simulating semiflexible inextensible fibers/polymers are based on discrete (bead-link or blob-link) models, which become prohibitively expensive in the slender limit when hydrodynamics is accounted for. In this paper, we develop a novel coarse-grained approach for simulating fluctuating slender filaments with hydrodynamic interactions. Our approach is tailored to relatively stiff fibers whose persistence length is comparable to or larger than their length and is based on three major contributions. First, we discretize the filament centerline using a coarse non-uniform Chebyshev grid, on which we formulate a discrete constrained Gibbs-Boltzmann (GB) equilibrium distribution and overdamped Langevin equation for the evolution of unit-length tangent vectors. Second, we define the hydrodynamic mobility at each point on the filament as an integral of the Rotne-Prager-Yamakawa kernel along the centerline and apply a spectrally accurate "slender-body" quadrature to accurately resolve the hydrodynamics. Third, we propose a novel midpoint temporal integrator, which can correctly capture the Ito drift terms that arise in the overdamped Langevin equation. For two separate examples, we verify that the equilibrium distribution for the Chebyshev grid is a good approximation of the blob-link one and that our temporal integrator for overdamped Langevin dynamics samples the equilibrium GB distribution for sufficiently small time step sizes. We also study the dynamics of relaxation of an initially straight filament and find that as few as 12 Chebyshev nodes provide a good approximation to the dynamics while allowing a time step size two orders of magnitude larger than a resolved blob-link simulation. We conclude by applying our approach to a suspension of cross-linked semiflexible fibers (neglecting hydrodynamic interactions between fibers), where we study how semiflexible fluctuations affect bundling dynamics. We find that semiflexible filaments bundle faster than rigid filaments even when the persistence length is large, but show that semiflexible bending fluctuations only further accelerate agglomeration when the persistence length and fiber length are of the same order.

Interplay between Brownian motion and cross-linking controls bundling dynamics in actin networks.

Morphology changes in cross-linked actin networks are important in cell motility, division, and cargo transport. Here, we study the transition from a weakly cross-linked network of actin filaments to a heavily cross-linked network of actin bundles through microscopic Brownian dynamics simulations. We show that this transition occurs in two stages: first, a composite bundle network of small and highly aligned bundles evolves from cross-linking of individual filaments and, second, small bundles coalesce into the clustered bundle state. We demonstrate that Brownian motion speeds up the first stage of this process at a faster rate than the second. We quantify the time to reach the composite bundle state and show that it strongly increases as the mesh size increases only when the concentration of cross-links is small and that it remains roughly constant if we decrease the relative ratio of cross-linkers as we increase the actin concentration. Finally, we examine the dependence of the bundling timescale on filament length, finding that shorter filaments bundle faster because they diffuse faster.

Simulations of dynamically cross-linked actin networks: Morphology, rheology, and hydrodynamic interactions.

Cross-linked actin networks are the primary component of the cell cytoskeleton and have been the subject of numerous experimental and modeling studies. While these studies have demonstrated that the networks are viscoelastic materials, evolving from elastic solids on short timescales to viscous fluids on long ones, questions remain about the duration of each asymptotic regime, the role of the surrounding fluid, and the behavior of the networks on intermediate timescales. Here we perform detailed simulations of passively cross-linked non-Brownian actin networks to quantify the principal timescales involved in the elastoviscous behavior, study the role of nonlocal hydrodynamic interactions, and parameterize continuum models from discrete stochastic simulations. To do this, we extend our recent computational framework for semiflexible filament suspensions, which is based on nonlocal slender body theory, to actin networks with dynamic cross linkers and finite filament lifetime. We introduce a model where the cross linkers are elastic springs with sticky ends stochastically binding to and unbinding from the elastic filaments, which randomly turn over at a characteristic rate. We show that, depending on the parameters, the network evolves to a steady state morphology that is either an isotropic actin mesh or a mesh with embedded actin bundles. For different degrees of bundling, we numerically apply small-amplitude oscillatory shear deformation to extract three timescales from networks of hundreds of filaments and cross linkers. We analyze the dependence of these timescales, which range from the order of hundredths of a second to the actin turnover time of several seconds, on the dynamic nature of the links, solvent viscosity, and filament bending stiffness. We show that the network is mostly elastic on the short time scale, with the elasticity coming mainly from the cross links, and viscous on the long time scale, with the effective viscosity originating primarily from stretching and breaking of the cross links. We show that the influence of nonlocal hydrodynamic interactions depends on the network morphology: for homogeneous meshworks, nonlocal hydrodynamics gives only a small correction to the viscous behavior, but for bundled networks it both hinders the formation of bundles and significantly lowers the resistance to shear once bundles are formed. We use our results to construct three-timescale generalized Maxwell models of the networks.

A fast spectral method for electrostatics in doubly periodic slit channels.

We develop a fast method for computing the electrostatic energy and forces for a collection of charges in doubly periodic slabs with jumps in the dielectric permittivity at the slab boundaries. Our method achieves spectral accuracy by using Ewald splitting to replace the original Poisson equation for nearly singular sources with a smooth far-field Poisson equation, combined with a localized near-field correction. Unlike existing spectral Ewald methods, which make use of the Fourier transform in the aperiodic direction, we recast the problem as a two-point boundary value problem in the aperiodic direction for each transverse Fourier mode for which exact analytic boundary conditions are available. We solve each of these boundary value problems using a fast, well-conditioned Chebyshev method. In the presence of dielectric jumps, combining Ewald splitting with the classical method of images results in smoothed charge distributions, which overlap the dielectric boundaries themselves. We show how to preserve the spectral accuracy in this case through the use of a harmonic correction, which involves solving a simple Laplace equation with smooth boundary data. We implement our method on graphical processing units and combine our doubly periodic Poisson solver with Brownian dynamics to study the equilibrium structure of double layers in binary electrolytes confined by dielectric boundaries. Consistent with prior studies, we find strong charge depletion near the interfaces due to repulsive interactions with image charges, which points to the need for incorporating polarization effects in understanding confined electrolytes, both theoretically and computationally.

Computational estimates of mechanical constraints on cell migration through the extracellular matrix.

Cell migration through a three-dimensional (3D) extracellular matrix (ECM) underlies important physiological phenomena and is based on a variety of mechanical strategies depending on the cell type and the properties of the ECM. By using computer simulations of the cell's mid-plane, we investigate two such migration mechanisms-'push-pull' (forming a finger-like protrusion, adhering to an ECM node, and pulling the cell body forward) and 'rear-squeezing' (pushing the cell body through the ECM by contracting the cell cortex and ECM at the cell rear). We present a computational model that accounts for both elastic deformation and forces of the ECM, an active cell cortex and nucleus, and for hydrodynamic forces and flow of the extracellular fluid, cytoplasm, and nucleoplasm. We find that relations between three mechanical parameters-the cortex's contractile force, nuclear elasticity, and ECM rigidity-determine the effectiveness of cell migration through the dense ECM. The cell can migrate persistently even if its cortical contraction cannot deform a near-rigid ECM, but then the contraction of the cortex has to be able to sufficiently deform the nucleus. The cell can also migrate even if it fails to deform a stiff nucleus, but then it has to be able to sufficiently deform the ECM. Simulation results show that nuclear stiffness limits the cell migration more than the ECM rigidity. Simulations show the rear-squeezing mechanism of motility results in more robust migration with larger cell displacements than those with the push-pull mechanism over a range of parameter values. Additionally, results show that the rear-squeezing mechanism is aided by hydrodynamics through a pressure gradient.

Langmuir-Scheaffer Technique as a Method for Controlled Alignment of 1D Materials.

A widely applicable method for aligning 1D materials, and in particular carbon nanotubes (CNTs), independent of their preparation would be very useful as the growth methods for these materials are substance-specific. Langmuir-Schaefer (LS) deposition could be such an approach for alignment, as it aligns a large number of 1D materials independently of the desired substrate. However, the mechanism and required conditions for alignment of 1D nanomaterials in a Langmuir trough are still unclear. Here we show, relying on numerical simulations of the Langmuir film compression, that the LS method is a powerful tool to achieve maximal alignment of 1D material in a controllable manner. In particular, 1D materials terminated with a suitable surfactant can align only if the velocity induced by the attraction between individual 1D entities is low enough relative to the flow speed. To validate this model, we achieved an efficient LS alignment of single-walled carbon nanotubes covered with a suitable surfactant relying on the numerical simulations. In situ polarized Raman microspectroscopy during the compression of Langmuir film revealed good quantitative agreement between the numerical simulations and the experiment. This suggests the applicability of the LS technique as a versatile method for the controlled alignment of 1D materials.

Zika virus dynamics: When does sexual transmission matter?

The Zika virus (ZIKV) has captured worldwide attention with the ongoing epidemic in South America and its link to severe birth defects, most notably microcephaly. ZIKV is spread to humans through a combination of vector and sexual transmission, but the relative contribution of these transmission routes to the overall epidemic remains largely unknown. Furthermore, a disparity in the reported number of infections between males and females has been observed. We develop a mathematical model that describes the transmission dynamics of ZIKV to determine the processes driving the observed epidemic patterns. Our model reveals a 4.8% contribution of sexual transmission to the basic reproductive number, R0. This contribution is too minor to independently sustain an outbreak but suggests that vector transmission is the main driver of the ongoing epidemic. We also find a minor, yet statistically significant, difference in the mean number of cases in males and females, both at the peak of the epidemic and at equilibrium. While this suggests an intrinsic disparity between males and females, the differences do not account for the vastly greater number of reported cases for females, indicative of a large reporting bias. In addition, we identify conditions under which sexual transmission may play a key role in sparking an epidemic, including temperate areas where ZIKV mosquito vectors are less prevalent.